Beyond single-photon localization at the edge of a Photonic Band Gap. 
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We study spontaneous emission in an atomic ladder system, with both transitions coupled near- 
resonantly to the edge of a photonic band gap continuum. The problem is solved through a recently 
developed technique and leads to the formation of a "two-photon+atom" bound state with fractional 
population trapping in both upper states. In the long-time limit, the atom can be found excited 
in a superposition of the upper states and a "direct" two-photon process coexists with the stepwise 
one. The sensitivity of the effect to the particular form of the density of states is also explored. 



I. INTRODUCTION 

The emergence of materials with photonic band gaps 
(PBG) Jl]j|, also referred to as photonic crystals, has 
stimulated a broad range of problems pertaining to the 
interaction of few-level atoms with unusual (structured) 
reservoirs. The unconventional photonic density of states 
(DOS) associated with such materials leads to a number 
of novel effects which pose severe demands on the theo- 
retical tools necessary for their prediction and/or inter- 
pretation. Certain basic issues related to such problems 
are encapsulated in the behavior of excited atomic states 
with transition frequencies around the edge of the DOS. 
Even the decay of an excited two-level atom (TLA) has 
been shown to exhibit novel behavior under such circum- 
stances, the most notable aspect of which is the forma- 
tion of the so-called "photon+atom" bound state (PABS) 
|@-p|. Even in the context of this simple arrangement, 
however, it has been difficult to handle the problem of 
more than one photon in the PBG reservoir. An ap- 
proach capable of surmounting this limitation has been 
presented recently in the context of a TLA and shown to 
provide the solution to the problem of the coupling of a 
TLA to both a PBG reservoir and a defect mode 
For a number of reasons, three-level atoms are of par- 
ticular interest in quantum optics and predictably their 
behavior in the context of structured reservoirs has been 
addressed. This includes ladder, lambda and V-type ar- 
rangements Mathematical difficulties, however, 
have limited the relevant investigations to the case of only 
one photon present in the structured reservoir. This is 
indeed rather limiting as to the scope of questions one is 
allowed to contemplate. 

It is our purpose in this paper to show how our ap- 
proach provides a way out of this limitation. Any of 
the above mentioned standard three-level systems can 
be addressed. We have chosen the present results on the 
ladder system which in open space involves a cascade of 



two photon emissions |l2| . As we show in the following 
sections, allowing both photons to be strongly coupled 
to the PBG reservoir we obtain what should be called 
"two-photon+atom" bound state. The formation of such 
a state is then found to be associated with a counterin- 
tuitive coherent evolution of the three atomic states. 
The first part of our study has been based on the so-called 
isotropic model for the DOS of the PBG reservoir Jl3| , 
which could be argued to effectively reduce the problem 
to one dimension in wavevector (fc)-space. A rather con- 
spicuous feature of that DOS is a sharp (divergent) peak 
at the edge of the gap. In order to explore the persis- 
tence of our predictions under more relaxed assumptions 
for the DOS, we have introduced a modified form rep- 
resented by a generalized Lorentzian profile. It is still 
isotropic, but it does not exhibit the singularity at the 
edge. The basic features of the predicted behaviour are 
still present with only some quantitative modifications. 
The paper is organized as follows: Section II contains 
a description of the system and a brief summary of 
the discretization method. In section III, we derive the 
equations of motion for the amplitudes involved in the 
wavefunction of the system Atom+Continuum, which are 
solved numerically. The dynamical behaviour of the sys- 
tem is thus investigated. In section IV the same problem 
is placed in the context of a generalized Lorentzian pro- 
file model for the DOS, while the results are summarized 
in section V. 

II. THE SYSTEM 

We consider a three level atom in a cascade con- 
figuration, with atomic levels |1),|2),|3) and energies 
hu)\, hu>2, respectively, where (u>i > U2) (Fig. |l]a). 
Both atomic transitions are coupled near-resonantly to 
the edge of a photonic band gap and are thus strongly 
modified. In the isotropic model and close to the band- 
edge, the dispersion relation of PBG materials, can be 
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approximated by the effective mass dispersion relation 
©0 



LU k ~ UJ e + A(k - fc ) 2 , 



(1) 



where u e , fco are the frequency and the wavenumber cor- 
responding to the band-edge and A is a material specific 
constant. The corresponding density of states (DOS) 
reads 



p(w) = 



Po 



:Q(UJ - LO e ), 



(2) 



where p is a material specific constant and O(x) the 
Heaviside step-function. 

Neglecting the zero-point energies of the field modes and 
adopting the rotating wave approximation, the Hamilto- 
nian of the system in the interaction picture (h = 1) is 
written as 

+^ 5 f(4a 32 e-^* - a x a 23 e iS ^), (3) 

A 

where a^i denote the atomic dyadic operators \k){l\ with 
k, le{l, 2, 3}; 6* = (u>i - lo 2 ) - w M , 8\ = (u 2 - w 3 ) - u\, 

while a^, a^{a\, a\) are the creation and annihilation op- 
erators of the structured continuum, which is coupled to 
the atomic transitions, via the respective coupling con- 
stants gp,g-p. 

The spectral response SR^^) corresponding to the PBG 
effective mass dispersion relation Eq. (Q) is given by 
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where the sum is over the polarizations, the integral runs 
over the solid angle and C\ is the effective coupling of the 
upper transition to the structured continuum . 
In order to deal with the double excitation into the struc- 
tured reservoir, we use the discretization method that we 
have presented recently ||. Briefly, we replace the den- 
sity of modes in Eq. (^|) near the atomic transitions (for 
lo < lj u ) by a collection of discrete harmonic oscillators, 
while the rest of the mode-density is treated perturba- 
tively since it is far from resonance. We propagate the 
total wavefunction (Atom+Modes) to obtain the dynam- 
ics of the system. 

For an arbitrary DOS, the frequencies of the discrete 
modes, are obtained through the relation 



uj 3+1 = Uj-x + 2/p(uj), 



(5) 



which for the case of the DOS of Eq. (g) can be reduced 
to 



ujj = uj e + j 2 Su>, (6) 

where Suj is chosen sufficiently small (Su) « 4.4 x 
10~ 4 C 2 ^ 3 ). The coupling C*' 1 ' of the upper transition to 
each one of the discrete modes is found using Eq. (^) 

(7) 
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where N is the number of discrete modes and lo u is the 
upper limit of the discretized part. A similar relation can 
be derived for the coupling of the lower transition, 
where now C2 will be the corresponding effective cou- 
pling. In general, the dipole moments (l|d|2) and (2|d|3) 
are different and thus for the rest of this paper we let the 
couplings C\ and C2 be different. 

Let us denote by |3, lj, l m ), a state of the combined sys- 
tem (atom+field) , where the atom is in state |3) and two 
photons have been emitted into the structured reservoir, 
populating the modes j and m respectively. Following 
this pattern of notation, the wavefunction of the com- 
plete system reads 



|*(*)) = o|l,0 1 0) 
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bj% l J ,0) + ^Q m |3,l J ,l m ). (8) 



III. EQUATIONS OF MOTION AND RESULTS 

The time dependence of the amplitudes, is governed 
by the Schrodinger equation and after eliminating the 
off- resonant modes [pfl, we find 
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where j, m are mode indices and for all discrete modes 
gW = Qf) = and g^ = gf = g£\ while the 
shift terms can be determined, converting the sum into 
an integral from to = u u to infinity. 

This set of equations is solved numerically for 150 dis- 
crete modes (u u w 10C 2 ^ 3 ), with the results presented 
in Figs. |||| We plot the population in the upper level 
(solid line), the intermediate level (dashed line) and the 
lower level (dot-dashed line), for various detunings of the 
transition frequencies from the band-edge. We define the 
relative detunings of the upper and lower transitions from 
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the band-edge as S12 — u>x—u>2~ l ^e and $23 = L03— lu c , 
respectively. 

In all figures, we can identify a "transient regime" , on a 

2/3 

short time scale of the order of C x , when part of the 
atomic population is lost. On a longer time scale ("dy- 
namic regime"), the populations in the atomic levels un- 
dergo oscillations, strongly dependent on the relative de- 
tunings from the band-edge, reflecting the emission and 
reabsorption of the photon(s). The localization at the 
atom of the photon emitted in the first transition, is ac- 
companied by emission of a photon in the lower transi- 
tion, which will be also localized at the atom. In analogy 
with the case of one photon in a TLA, two photons are 
now backscattered to the atom after tunneling a charac- 
teristic distance and reexcite it. 

As is known from the coupling of a TLA to the PBG 
reservoir, the dressing of the atom by its own radiation 
causes splitting of the atomic levels. This splitting is suf- 
ficiently strong to push one level of the doublets outside 
the gap and the other inside. The dressed state outside 
the gap looses all its population in the long-time limit, 
while the one inside the gap is protected from dissipation 
and thus is stable. The number of stable localized states, 
is intimately connected to the behavior of the system in 
the long-time limit. One such state gives rise to steady 
state population in the excited level, while two of them 
lead to an undamped beating of the system between the 
two non-decaying states 

For our three-level atom, two transitions are at the 
edge, and thus more than one stable localized state 
can be found in the gap. In analogy to the "single- 
photon+atom" bound state, we have the formation of 
a "two-photon+atom" bound state, which exhibits pop- 
ulation trapping in both excited states, in the long-time 
limit. Thus, the atom is finally excited, in a superposition 
of the upper states (Fig. ||). This is a novel behavior, due 
to the fact that both transitions and not only one, are 
coupled to the same structured continuum. Note that, 
even in the case that only the (|2) — * |3)) transition is 
at the edge of the gap, it is the intermediate level that 
exhibits non-zero steady-state population but not the up- 
per level §0. 

In the language of dressed states, the oscillations in the 
populations of the atomic levels reflect the interference 
between the dressed states of the atom. We addition- 
ally note that the oscillations in the populations of the 
intermediate and lower level are "in phase". This is a 
rather surprising result, since one would expect the pop- 
ulation in the intermediate state to be maximum when 
the population in the lower state is minimum and vise- 
versa. The phenomenon is even more pronounced in Figs. 
^. After an initial transient regime in which about 15% 
of the population is lost, the remaining population oscil- 
lates between the upper level and the two lower levels, in 
a non-dissipative way. The undamped oscillations imply 



the beating of the system, between more than two non- 
decaying dressed states. 

The "in phase" oscillations, stem from the localization 
of both photons at the site of the atom. The system 
can make the transitions |1) <-» |3) either with a stepwise 
process (|1) <-* \2) *-* |3)) or with a two-photon process 
(|1) |3)). Which of the two routes the system will fol- 
low to arrive at |3) depends on the detuning 612- 
Specifically, if the upper transition is outside the gap, 
whatever the detuning from the edge of the lower transi- 
tion is, it seems preferable for the system to "decay" via 
the stepwise process rather than the two-photon process 
indicated in Fig. 4. This is no different from the behaviour 
of a ladder system in open space. On the contrary, if the 
upper transition is inside the gap, the system evolves in 
time as if the upper state were coupled to the intermedi- 
ate state via a single-photon process, and simultaneously 
to the ground state via a "direct" two-photon process, 
with respective frequencies fii and f^. This can not be 
anticipated on the basis of the behaviour in open space 
and it is what we meant by counter-intuitive in the in- 
troduction. 

To gain further insight into this effect, we can adopt 
a simple 3-level model without dissipation and assign 
a single-photon Rabi frequency fl\ between |1) and \2) 
and a two-photon Rabi frequency O2 between |1) and |3) 
(Fig. |l|b). An analysis through standard rate equations 
shows that the populations of the atomic levels oscillate 
with the same frequency £1 = y'fq + fif, with the oscil- 
lation of the two lower levels being in phase. The ratio 
of the amplitude of the oscillation of the intermediate 
level to that of the lower level is related to the ratio of 
the Rabi frequencies, i.e. Sli/0,2- In Fig. ^|b represent- 
ing the result of the numerical calculation for the system 
in the PBG reservoir, we note that the dashed and the 
dot-dashed lines, corresponding to the populations in the 
intermediate and ground levels, respectively, are indistin- 
guishable. This implies that the corresponding effective 
Rabi frequencies f2i and £^2 are practically equal. It is 
the combination of the coupling constants and detunings 
that conspire to produce that behaviour. 
The effective detuning of the |3) «-> |1) transitions from 
the band-edge for the "direct" two-photon process is de- 
fined as A^ 2 ) = cl>i — CJ3 — 2lu e . From the known dynamics 
of a TLA with transition frequency at the edge of the gap, 
depending on A^ 2 \ we may expect suppression or inhi- 
bition of the "direct" two-photon process. Specifically, 
for detunings inside the gap (A^ 2 ) < 0), the "direct" 
two-photon emission should be totally or partially sup- 
pressed, while for detunings outside the gap and almost 
at the edge (A^ 2 ) > 0), it should be enhanced due to the 
high density of final available states. 
For transitions symmetrically placed around the band- 
edge (5 12 — —S23), with the upper one being inside the 
gap (6x2 < 0), the detuning for the "direct" two-photon 
transition is exactly at the edge (A^ 2 ' = 0) (Fig. ||). 
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For the parameters used in Figs. || and without tak- 
ing into account the "direct" two-photon process, the 
"single-photon+atom" bound state for the upper state 
should be metastable, in the sense that the main part of 
the population should be lost in the long-time limit. On 
the contrary, we find that this does not happen and the 
part of the population that has not been lost in the tran- 
sient regime, oscillates between the upper and the ground 
state. These oscillations are not reflected in the interme- 
diate level's population which remains almost constant 
with some oscillations of negligible amplitude. This be- 
havior definitely indicates the coupling of the upper level 
to the ground level via a "direct" two-photon process as 
described above. Note again the "in phase" oscillations 
for the two lower levels. 

Choosing S 23 such that A^l > (Fig. |), the mam 
part of the population is indeed lost in the long-time 
limit. The difference between Fig. [5|a and Fig. [6] is 
the detuning of the lower transition from the band-edge. 



In both cases (623 = 2C. 
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the behav- 



ior of a TLA in the long-time limit is the same i.e. the 
population is lost. For the Ladder system, however, we 
note an oscillatory behavior where part of the popula- 
tion is trapped to the atom in the long-time limit for 



(523 = 2C^ 3 and complete decay for 823 — AC^' - The 
photon-atom bound state formed due to the upper tran- 
sition becomes therefore metastable as soon as S23 is 
chosen such that A' 2 ) > 0. This is a conclusion that has 
been checked for various detunings of the atomic transi- 
tions from the band-edge but with the upper one always 
in the gap (612 < 0). 
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IV. GAP WITH A LORENTZIAN PROFILE OF 
DOS 

The isotropic model for the DOS of a photonic crys- 
tal has been employed quite extensively in the literature. 
That is the model we have employed in the previous sec- 
tions in order to have a direct assessment of the new 
results vis a vis those obtained in previous work with 
only one photon in the reservoir, where the same DOS 
has been assumed. It is, however, well known that this 
DOS is somewhat artificial, with a divergence as the fre- 
quency approaches the band edge. The infinitely high 
sharp peak at the edge tends to exaggerate many effects 
and it is important to check predictions made on the ba- 
sis of this model against other forms of DOS. With this 
intention, we introduce here a DOS which although still 
isotropic, does not exhibit a divergence at the edge. 
What is essential for an appropriate model of the DOS is 
that it exhibit a dip and also that it tend to the open 
space DOS as the frequency becomes much larger or 
smaller than the mid-gap frequency. We have chosen to 
adopt as a model of such a DOS an inverted Lorentzian 



of higher order given by the expression 



PM = Pa 
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pr 



(cj - Lu ) n + r« 



(13) 



First of all, this DOS approaches the open space value 
Po for \w — ljJ o\ >> T. Second, it does not exhibit a di- 
vergence at the edge. In fact, the "edge" is not infinitely 
steep, but does rise more steeply, as n — > 00. 
It could be argued that this DOS does not exhibit a clear 
edge and possesses a zero only at one point. It should be 
kept in mind on the other hand that, in a realistic PBG 
material, the gap does not necessarily mean a true zero 
but a range of frequencies over which the DOS is several 
orders of magnitude smaller than that of open space [^5) . 
Taking n sufficiently large, in Eq. (|ll|), one can obtain a 
range of frequencies over which p(lo)/pq is smaller than 
a desired value. For n = 6, for example, p(uj)/po < 10~ 6 
for lu — lu ± O.ir. One can further combine the inverted 
Lorentzian with step functions in order to simulate a true 
zero over a range of oj, if so desired. The shape of the 
DOS given by Eq. (|l3|), can be viewed as a compromise 
between the isotropic and the anisotropic models [13|, 
and has been used in the literature for n = 2 |l^ , [l6| , 17|. 
It should be stressed once more that the modified DOS 
in no way relaxes its isotropic in k space nature. It only 
eliminates the singularity. 

We proceed now to the exploration of the ladder system 
under the DOS of Eq. ( |l3| ) adopting the specific case 
of n = 6. It is one of the strengths of the technique 
of discretization that it can be implemented with essen- 
tially any DOS, provided the manner of discretization is 
adapted to the demands of the particular form. For the 
case under consideration, which is symmetric around lu , 
we chose a range loi ow < uj < uj up within which the DOS 
is replaced by a sufficiently large number (in this case 
150) equidistant discrete modes. The spectral response 
corresponding to Eq. (|13]) for n — 6 is of the form 



SR-M = I 
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(14) 



where 71 is the decay rate of the upper state in free space. 
The coupling of the upper transition to the i th discrete 
mode is frequency dependent and is given by 




r 6 



Au, (15) 



where Auj ps O.2772 is the spacing between two discrete 
modes and uii is the frequency corresponding to the i th 
mode. An analogous relation gives the coupling of the 
lower transition, Q^ 2 \ to the i th mode, where 71 is re- 
placed by 72, the free space decay-rate of the interme- 
diate state. By way of comparison illustrating the phi- 
losophy of the discretization procedure, note that here 
the mode-spacing is uniform but the coupling constant 
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frequency-dependent; while in section II (Eqs. g and [7]) 
it is the other way around. 

The rest of the mode-density for u> > u> up and to < uji ow , 
is treated perturbatively M leading to a shift for the 
two upper levels. These shift terms differ from those 
in Eqs. (p[-[L2|) and are not given explicitly. They can, 
however, be determined numerically. The relative posi- 
tions of the upper and lower transition frequencies are 
now defined with respect to the central frequency lu of 
the gap, in terms of the detunings: S12 = oj\ — L02 — lu 
and S23 — UJ2 — — oj , respectively. In Fig. [?], we 
present the evolution of the population in the states of 
the Ladder system as a function of time, for a particular 
combination of detunings. As in the previous figures cor- 
responding to the isotropic model, both upper and lower 
levels exhibit non-zero steady-state population, as a con- 
sequence of the "two-photon+atom" bound state. The 
oscillations, however, in the atomic populations, which 
can be interpreted as interference between the dressed 
states, are not present. It seems that this oscillatory be- 
havior is strongly related to the isotropic model. For the 
Lorentzian profile, the part of the doublet that is pushed 
outside the gap, decays much faster than the isotropic 
model would predict. From this we conclude that in the 
isotropic model, the "dynamic regime" which follows the 
initial transient regime and is dominated by the emis- 
sion and reabsorption of photons, or else the interference 
between the various dressed states, is much more pro- 
nounced than in the Lorentzian model, or we would ar- 
gue any model that does not exhibit the highly peaked 
feature of the isotropic model . Without the concen- 
tration of the density of states around a peak, the notion 
of dressed states is diluted. Nevertheless, the coherent su- 
perposition of states and the "two-photon+atom" bound 
state persists even in this model. Atomic populations re- 
main trapped for long times, long compared to 1/72- 
A further issue needs to be brought up here, in connec- 
tion with the time scale of the persistence of any "pho- 
ton+atom" bound state. As long as the model for the 
DOS involves an exact zero over some frequency range, 
the "photon+atom" bound state lives for ever, which 
mathematically implies a non-zero fractional population 
trapping in the limit t — > 00. In reality, however, the 
DOS, more often than not, will not involve an exact zero 
but a deep minimum, as already mentioned in the be- 
ginning of this section. As a result, the life-time of the 
"photon+atom" bound state, will be long on some time 
scale but not necessarily for t — > 00 in the mathematical 
sense. It is this finiteness of the lifetime of the "two- 
photon+atom" bound state that would cause the popu- 
lations of states |1) and |2) in Fig. 7 to decay to zero, if the 
calculation were extended to sufficiently long times. The 
essential point therefore is that the population may re- 
main trapped for such a long time which, for all practical 
purposes, is equivalent to t — > 00. 



V. SUMMARY 

We have investigated the dynamics of a Ladder atomic 
system with both transitions coupled to the same struc- 
tured reservoir. This has been possible through a dis- 
cretization approach according to which the density of 
modes over a range near the atomic transitions is re- 
placed by an appropriately chosen collection of discrete 
harmonic oscillators, with the rest of the mode-density 
treated perturbatively. We have found that this system 
supports a "two-photon+atom" bound state which leads 
to a fractional population trapping in both of the upper 
states and the atom can be in a superposition of the up- 
per levels even in the long-time limit. In the presence of 
the two photons at the site of the atom, we have shown 
that the atom has two paths for the |1) |3) transition, 
and found that a "direct" two-photon process coexists 
with a stepwise one. Which of the two dominates is de- 
termined mainly by the detuning of the upper transition 
from the band-edge. We have further explored the persis- 
tence of this effect under much more relaxed forms of the 
DOS and shown that, although quantitatively modified, 
the basic effect remains. 
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(a) (b) 

FIG. 1. (a) Schematic representation of the atomic system 
and the possible transitions, (b) The upper state simultane- 
ously coupled to the intermediate state and the ground state 
via a single- and two-photon process respectively. 




0.0 10.0 20.0 30.0 40.0 

Time 



FIG. 2. The population in the atomic states as function of 
time. The solid line is for |1), the dashed line for |2) and the 
dot-dashed line for |3). Parameters: C2 = 1.5&, £12 = — C^J Z 
and &2Z =0. The time is in units of 
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FIG. 3. The population in the atomic states as function 
of time. The solid line is for |1), the dashed line for 2) 
and the dot-dashed line for |3). Parameters: C2 = 1.5&, 
(a) <5i2 = -2Ca /3 and 823 = 1C 2/S ; (b) 612 = -2C 2 J 3 and 
<5 2 3 = 1C'^ 3 ' ■ The time is in units of C^ 3 . 



FIG. 4. The population in the atomic states as function 
of time. The solid line is for the dashed line for 2) 
and the dot-dashed line for |3). Parameters: C2 = 1.5&, 
(a) S 12 = 1C 2 2 ' 3 and 823 = -IC 2 / 3 ; (b) 8 12 = 2C 2 2 ' 3 and 
$23 = iC 2 ^ 3 . The time is in units of C 2 ^ 3 . 
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FIG. 6. The population in the atomic states as function of 
time. The solid line is for |1), the dashed line for \2) and the 
dot-dashed line for |3). Parameters: C 2 = 1.5&, 612 = —2C^ 3 
and 623 = 4C 2 2/3 . The time is in units of C 2/3 . 



FIG. 5. The population in the atomic states as function 
of time. The solid line is for |1), the dashed line for 2) 
and the dot-dashed line for |3). Parameters: C2 = 1.5&, 
(a) <5i2 = -2C 2 / 3 and S 23 = 2C^ Z ; (b) 612 = -2C\ ! ' A and 
&23 = 2C 2 ^ 3 . The time is in units of C 2 ^ 3 . 
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FIG. 7. The population in the atomic states as function 
of time. The solid line is for (1), the dashed line for |2) and 
the dot-dashed line for |3). Parameters: 71 = O.572, T = 72, 
S12 = O.I72, S23 = O.372 and u>up = —uiow = 2O72. The time 
is in units of 72 . 
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